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The study of systems of strongly correlated electrons has received an increasing 
attention in the last years, especially after the discoveries of the Fractional Quantum 
Hall Effect and the High Temperature Superconductivity. In this context, one of the 
most studied is the Hubbard model [1] which describes the dynamics of non-relativistic 
electrons moving on a lattice. Its Hamiltonian contains a kinetic term and a quartic 
potential term which effectively describes the Coulomb interactions among electrons. 
Another lattice model that has been suggested as an appropriate starting point 
■ for the theory of the High Temperature Superconductivity [2,3], is the so-called t-J 
model [4] which is characterized by the absence of doubly occupied electron states. 
Its Hamiltonian contains a kinetic term of strength t responsible for the hopping 
of electrons from one lattice site to its nearest neighbors, and a potential term of 
strength J which describes nearest neighbor spin exchange interactions. If J « t, 
, the t-J model can be mapped onto a Hubbard model with very strong Coulomb 
repulsion; for other values of the coupling constants instead, the two models exhibit 
a quite different behavior. A particularly interesting feature of the t-J model is the 
J> ■ presence of a supersymmetry [5]: In fact, when J = ±2t, the Hamiltonian of the t-J 
model commutes with the generators of the superalgebra 5'L'^(1|2) [6]. It is remarkable 
that precisely at the supersymmetric point J = ±2 1, the t-J model in one dimension 
is exactly solvable by Bethe Ansatz [7,8]. We recall that the supersymmetric t-J 
^T) [ model cannot be derived from the Hubbard model, but instead can be mapped onto 
^ ■ a quantum lattice gas of hard core bosons and fermions. This latter system was 
formulated in the mid-seventies and solved exactly by Bethe Ansatz in one dimension 

D , The availability of a complete solution for the one-dimensional supersymmetric 

^ ' t-J model allows a detailed analysis of its collective excitations [10]. In particular two 
kinds of low lying quasiparticle excitations have been identified above the antiferro- 
^ ■ magnetic ground state (J > 0): One carries charge and no spin and is called holon; the 
other carries spin and no charge and is called spinon. This remarkable phenomenon 
of spin-charge separation, which is rigorously proven only in one dimension, has been 
conjectured to occur also in two dimensions [3,11]. To implement this separation be- 
tween spin and charge, one often uses the formalism of the slave operators [12] and 
makes the Ansatz in which the electron is a product of a charged spinless antiholon 
and a neutral spinon of spin 1/2. Since the electron is a fermion, holons and spinous 
must have complementary statistics. For example, the holon may be fermionic and 
the spinon bosonic (slave fermion representation [13]), or vice- versa the holon may be 
bosonic and the spinon fermionic (slave boson representation [3,11]). Actually, since 
in two dimensions the statistics is arbitrary [14,15], there are in principle other possi- 
bilities; in particular it has been suggested [16] that holons and spinous are semions, 
i.e. anyons that are exactly half-way between bosons and fermions. In [17,18] it has 
been shown that these different statistics assignments can be interpreted in a unified 
way within the framework of the bosonization procedure of two-dimensional systems. 
More precisely, the slave boson and slave fermion representations arise as two different 
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gauge fixings in the abefian bosonization of the t-J model, whereas the semion case 
naturally occurs in the non-abelian bosonization. 

In this paper we present a generalization of the analysis of [17,18], and in partic- 
ular we show that performing an abelian bosonization with two independent Chern- 
Simons fields, the holons and the spinous of the two-dimensional t-J model may 
become anyons of arbitrary complementary statistics, in such a way that the product 
of one antiholon and one spinon (i.e. the electron) is a fermion. We call this the slave 
anyon representation, which interpolates continuously between the slave boson and 
the slave fermion representations; in this way the semionic statistics of [16] simply 
becomes a particular case. Therefore, we deduce that at the formal kinematical level, 
all consistent statistics assignments are possible and only a dynamical analysis can 
eventually select one among these many possibilities. 

Another interesting fact is that the supersymmetry of the t-J model at J = ±2t 
becomes simpler and more transparent when formulated in terms of holon and spinon 
operators. In fact, since these have complementary statistics, it is natural to expect 
that the fermionic supercharges transform holons into spinous and vice- versa ^ . Indeed 
this is what happens, because the superalgebra SU{1\2) is linearly realized on holons 
and spinous in contrast to the non-linear realization in terms of the electron operators. 
In this paper we discuss in particular the supersymmetry of the t-J model in the slave 
anyon representation and analyze the braiding properties of the supercharge densities 
with holons and spinous. 

The paper is organized as follows: In Section 2 we review the main features of 
the t-J model, write its Hamiltonian and present its bosonization with two indepen- 
dent Chern-Simons fields to yield holons and spinous with complementary statistics. 
In Section 3 we discuss the supersymmetry of the t-J model in the slave fermion 
formalism, whereas in Section 4 we introduce slave holons and spinous in the anyon 
representation and discuss their braiding properties. Then we prove the invariance un- 
der SU{1\2) of the Hamiltonian of the t-J model at the supersymmetric point J = 2t 
in the slave anyon representation. Finally, Section 5 presents our conclusions, and 
Appendix A contains an alternative proof of the supersymmetry of the Hamiltonian 
of the t-J model when J = 2t. 



2. The t-J Model and its Bosonizations 

The (two-dimensional) t-J model describes a lattice system of strongly correlated 
non-relativistic electrons. It was proposed a few years ago by Zhang and Rice [4] as 



^ To avoid confusions, we stress that here the terms holons and spinous are used 
in a different way with respect to [10] where holons and spinous denote different 
supermultiplets of SU {1\2). 
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an effective model for the electronic dynamics in the copper-oxide planes of the high 
temperature superconductors. Given a two-dimensional square lattice O of spacing 
a = 1, one assumes that in each site at most one electron can be accomodated, i.e. 
doubly occupied electron states are excluded from the physical Hilbert space of the 
system; this is a way to simulate strong Coulomb repulsion. The t-J model with 
chemical potential fi is described by the following Hamiltonian 



Htj = V 



E E (4Wca(j)+h.c.) 

<i,j> a=±l 

+J J2 (s(^)■S{j)-\n{^)n{j)]+^^J2n{^) 

<i,j> ^ ^ ien 



(2.1) 



where ^ denotes the sum over nearest- neighbor lattice sites; cl^{i) and Ca{i) are 

<i,J> 

respectively creation and annihilation operators for electrons with third component 
of the spin q;/2 at the site i. They satisfy the usual anticommutation relations 

{ca(i),4{j)} = 6^^S{i,j) (2.2) 

where S{i,j) is the lattice delta function. The vector S(i) is the spin operator whose 
components are 

S'-ii)^ E '^^'^ a = 1,2, 3, (2.3a) 

((j" being the Pauli matrices), and n{i) is the number operator 

l(i)= J2 "a(^)= E 4(^)Ca(i) . (2.3b) 



a=±l a=±l 



Finally, V denotes the Gutzwiller projection operator which ensures that doubly oc- 
cupied electrons states are removed from the Hilbert space of the system [19]. The 
parameters t and J in (2.1) are coupling constants which we take to be positive: 
t denotes the strength of the kinetic nearest-neighbor hopping term; J denotes the 
strength of the spin exchange nearest-neighbor interaction. We recall that for J « t 
the hamiltonian (2.1) can be derived from the Hubbard model in the limit of (infi- 
nite) strong on-site repulsion and that for t = it reduces to the antiferromagnetic 
Heisenberg XXX model [1]. 

In the following we will discuss the bosonization procedure with Chern-Simons 
fields and the presence of anyons, which are specific features of the two-dimensional 
t-J model; however, most of the considerations that we present hereinafter are valid 
for lattices of any dimensions. As we mentioned, doubly occupied states are excluded 
from the Hilbert space of the t-J model; therefore if we denote by |0) the Fock vacuum 
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which satisfies Cq, (z)|0) = for all i E O, and a = ±1, there are only three possible 
electronic states at any given lattice site i, namely 



|0) , clim , cl(z)|0) . 



(2.4) 



The fourth state c^(?') c'_(?')|0), which would be allowed by Fermi statistics, is removed 
by the Gutzwiller projector V. There is an equivalent way to realize this truncation 
in the Fock space: Instead of using the operators cj^, to construct the physical states 
and then the Gutzwiller projector to remove those with double occupancy, one can 
use the projected fermionic operators 



c'a{i) = (l - ^-a(^)) Ca{i) , 

c'i{i)^cUi) (l-n_«(z)) 



(2.5) 



for a — ±1. It is immediate to verify that these operator satisfy a hard-core constraint 
even when their spin index is different. Indeed one has 



c'^ii) c'f,{i)=c'l{i) 4{i) = 
for all a and Furthermore, using (2.2) it is straightforward to verify that 

{c'a{i),c'^{j)} = {c'i{i),4{j)]=0 , 

{c'aii) , 4{j)} = 5(^, j) (l - \n{i)\ + X^5"(z)< 

^ ^ a=l 



(2.6) 



(2.7) 



The operators S°'{i) and n{i) are defined respectively in (2.3a) and (2.3b), but it is 
easy to realize that they can be written indifferently in terms of Ca or c' a- In fact, 
with simple algebra we see that 



a/3 



c'/3(^)= E 



a,/3=±l 



cp{i) = S''{i) , (2.8a) 



a/3 



and 



2 

n'a{i) = 4ii) c'aii) = ci{i) Ca(*) (l - n^aii)^ = riaii) , 



(2.8b) 



where the last equality follows from na{i) ri-aii) = which is a consequence of the ab- 
sence of doubly occupied states. Using the projected operators (2.5), the Hamiltonian 
(2.1) becomes 



tJ 



E E (c'iwcaOO+h-c.) 

<i,j> a=±l 

+ J E (s(z)-s(j)-^^w^O'))+/^E^(^) 



<^,J> 



(2.9) 



Since the Hilbert space of the t-J model does not contain any doubly occupied 
states, the hopping term of the Hamiltonian can be dropped at half-filling, i. e. when 
the number of electrons is equal to the number of lattice sites. In this case the t-J 
model becomes equivalent to the (antiferromagnetic) XXX model. The situation is 
clearly different below half-filling: In such a case there are some empty points and 
some electrons can hop into them from neighboring occupied sites with a rate that is 
proportional to the coupling constant t. A deviation from half-filling is usually called 
"doping" and can be realized by suitably choosing the value of the chemical potential 
H in the Hamiltonian. Several authors suggested that, at small doping, the spectrum 
of the lowest lying excitations above the antiferromagnetic ground state consists of 
holons and spinous [2,3,11]. The holons carry charge and no spin, whilst the spinous 
carry spin and no charge. Only an accurate study of the dynamics could allow to 
check if this mechanism of spin-charge separation actually occurs. To this aim it may 
be useful to perform a change of variables and write the projected electron operators, 
which obviously have both charge and spin, as products of an antiholon and a spinon. 
More precisely one writes 

c'a(i) = e\i) s^{i) , c'i(i) = e{i) 4(i) (2.10) 

where e^(z) and e(z) are respectively the creation and annihilation operators for 
charged spinless holons, and s]^{i) and Sa(i) are respectively the creation and an- 
nihilation operators for neutral spinous of spin a, and are such that 

et(^)|0) = 5«(^)|0) = 

where |0) is the Fock vacuum of Eq. (2.4) representing a hole. The requirement that 
there are no doubly occupied electron states, translates into the constraint 

4W^a(^)-l = (2.11) 

a=±\ 

for all z G O. This constraint implies that in each lattice site one of the following 
three possibilities must occur: Either there is one hole or there is one electron with 
spin -1-1 or one electron with spin —1. There are no other possibilities, just like in the 
original electronic description where there were only three allowed configurations at 
each site. 

Another important issue is the statistics of holons and spinous. Since the elec- 
tron operators must be fermionic, e{i) and ««(«) cannot have the same statistics. In 
principle we have two possibilities: The holon is a fermion and the spinon is a boson, 
or vice-versa, the holon is a boson and the spinon is a fermion. These two scenarios 
are usually called the slave-fermion and the slave-boson descriptions respectively, and 
both of them have been proposed and studied in the literature (the slave-fermion 
approach for the t-J model has been considered in [13,20], while the slave-boson ap- 
proach has been considered in [3,11,21]). However, for two-dimensional systems there 
are other possibilities, since in two dimensions the statistics is not necessarily bosonic 
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or fermionic but can be anyonic [22] (for a review see for instance [14,15]). Indeed, in 
[16] Laughlin conjectured that the two-dimensional holons and spinous are semions, 
i.e. anyons that are exactly half-way between bosons and fermions. Actually, as 
we will show later, holons and spinous for the t-J model in two dimensions can in 
principle be anyons of arbitrary statistics with the only constraint that the composite 
object formed by one holon and one spinon be a fermion. 

A few remarks are in order at this point. First of all, we should keep in mind 
that the replacement of the electron operators with products of holons and spinous 
according to (2.10) is a purely formal operation. Differently from the electron op- 
erators c^(i) and c2{i), the holon and spinon oscillators do not commute with the 
constraint operator </)(z); therefore the spin-charge separation, i.e. the excitation of 
a single holon or spinon, can occur only if the constraint (2.11) is somehow released. 
Secondly, the issue of the statistics of holons and spinous cannot be resolved at this 
formal stage, since all the above mentioned assignments of statistics are formally cor- 
rect. In fact, to make any progress it would be necessary to determine what is the 
dynamics of holons and spinous and to find their Hamiltonian. The first step towards 
this aim is simply to substitute (2.10) into the Hamiltonian of the t-J model. For 
definiteness we choose the slave fermion description, but similar considerations would 
apply equally well to the slave boson formalism. In the slave fermion description, the 
fermionic holon operators obey 



{eii),eHj)} = Sii,j) , 
and the bosonic spinon operators obey 



(2.12a) 



(2.12b) 



with all other possible (anti) commutators vanishing. A straightforward substitution 
of (2.10) into the Hamiltonian (2.9) leads to ^ 

iytj=| + t E {eHj)e{i)si(i)sM + ^-c.)+j E eWen^)e(j)et(j) x 

I <i,j> <i,j> 

' 3 

i^^ocii) <^ Spii) 4'(i) <'/3' Sp'ij)) - 4(«) Sa{i) sl,{j) Sa'U) 

.0=1 

L e(i) et(i) sj,(i) sa(i) + i A(i) 0(i) 



(2.13) 

where we have introduced the Lagrange multiplier A(z) to enforce the constraint 
(2.11). 



^ To simplify the notation, here and in the following we understand the summation 
symbols over repeated spin indices. 
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Using the anticommutation relations of the holons, the J-term of (2.13) can be 
rewritten as follows 



J 

4 



<i,j> 



[l-et(i)e(i)-et(j)e(j) + et(i)e(i)et(j)e(j) x 

3 

(4W<^S/3Wsi'(iX'/3'S/3'0')) - si{i)Saii)s'l,{j)Sa'ij) 

' (2.14) 

Due to the constraint (2.11), it is easy to realize that only the 1 of the first square 
bracket of (2.14) yields a non vanishing contribution. In fact the other three terms, 
when multiplied by the second square bracket, give rise to expressions containing 
annihilation operators in the same point both for holons and spinous, and hence 
vanish when applied to physical states satisfying the constraint (2.11). With a similar 
argument, we can also rewrite the chemical potential term simply as 

Thus, the Hamiltonian (2.13) can be drastically simplified and one gets 
Hu=J2 \+^ (e^(i)e(z)4(0«a(j)+h.c.) 

<i,j> I 

J ^ 

+ I </3 '^^^ ^«'(-^') <'I3' (•?■)) 

a=l 

-si{l) Sa{i) sl,{j)Sa'{j)\ 

(i)+iA(i)0(i) . 

We stress that the formal manipulations that yielded the Hamiltonian (2.15) are valid 
in any dimensions. In the particular case of two dimensional lattices, it is possible to 
obtain (2.15) in an alternative way based on the abelian bosonization, as shown in 
[17,18]. The two-dimensional abelian bosonization consists of replacing the original 
fermionic fields by new bosonic degrees of freedom interacting with an abelian Chern- 
Simons gauge field [23] . In this framework, holons and spinous arise, roughly speaking, 
as the modulus and phase of the field which bosonizes the electron operators, and thus 
at first are both bosonic . The final statistics assignments are dictated by the choice 
of the coupling constants between the bosonic holons and spinous on the one hand, 
and the abelian Chern-Simons field on the other. 



(2.15) 



^ Actually the correspondence is more involved than this, and moreover one of 
the two fields must be a hard-core boson. For details we refer to the original papers 



[17,18]. 
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We now give some details of this construction following [18], and later we present 
its generalization. In the formalism of second quantization, the gran canonical parti- 
tion function for the t-J model at temperature /3 and chemical potential ^ is written as 
an Euclidean functional integral over a complex Grassmann field tl^a with antiperiodic 
boundary conditions, according to 



Q(/3,/x) = j ViP e-^W 



where 



dr 



(2.16) 



(2.17) 



In this formula Htj{ip) denotes the functional obtained from the Hamiltonian (2.9) by 

„ r x / /.N „„J „/t /.N „.:.i- 



111 imti iormuia riuyipj uenoies lue lunciioncii ouiameu irom ine namniommi yz^.i 
replacing the fermionic operators c' oi{i) and c')j,(i) with the Grassmann fields V'a(^, t) 
and 'ip^{i,T) respectively. This functional is explicitly given by 

<i,j> 

where 



(2.18a) 



3 

(^«(^' ^) <I3 ^) C'O', r) al,p, (j, r)) 

,o=l 

- (i, r) V'a (i, r) V'a' (i, ^) V^a' (j, ^) 



(2.18b) 



As is well known, the path integral over anticommuting Grassmann fields automat- 
ically takes into account the sign factors due the fermionic statistics. An equivalent 
expression for Q(/3, can be obtained by using, in place of •i/'a(^, t) and "fpaii, r), com- 
plex commuting hard core bosonic fields ipa{i-, ^) and (/'^(i, r) together with a suitable 
prescription to implement the original fermionic statistics. For two dimensional sys- 
tems, this prescription is quite simple and elegant: One couples these bosonic fields 
to an abelian Chern-Simons gauge field A\ (A = 0, 1, 2) with level k = 1, and averages 
over all its possible configurations [24]. More precisely, the following identity holds 



Q(/3,//) = V^j e-^(^) = ( Vipe 



(2.19) 



fe=i 



where S{(fi; A) is the bosonized action obtained from S{i(^) by replacing i/j^ with (p^i 
minimally coupled to A\. The symbol ( )k denotes the expectation value over the 
Chern-Simons field, namely 



0{A)), 



JVA e-^cs(^) 



(2.20) 
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where 0{A) is any functional of A and 



Shsi^) = ^jyTj d^x e^'^MA(T, x) d^A^ir, x) (2.21) 

is the abelian Chern-Simons action of level k. 

The bosonization identity (2.19) is essentially based on the fundamental property 
that the expectation values of Wilson loops in Chern-Simons theories are topological 
invariants [24,25]. To see the implications of this fact, let us consider an arbitrary 
braid b of N objects corresponding to the evolution of particles between their 
configurations at Euclidean times r = and t = /3. Due the periodic boundary 
conditions on r, the braid b actually forms a closed link L{b), and an abelian Wilson 
loop can be defined as follows 

2 

w(^L{b),A)=e (2.22) 

where the line integral in the exponent is computed alond the link L(6) with a suitable 
choice of framing [25]. If A\ is a Chern-Simons gauge field with action (2.21), one 
can prove that 



where 



(^L(6),An =e^^"(^) (2.23) 
/ k 



n(6) = ^[n+(6)-n_(6)] (2.24) 

is the the winding number of the braid 6, i.e. the difference between the number of 
overcrossings n-|-(6) and the number of undercrossings n_(6). If /c = 1, we have 

. ^ V.=i (2.25) 

= giTT (n+(6)-Fn_(6)) ^ 

where cT(7r(6)) is the signature of the permutation 7r(6) associated to the braid b. From 
(2.25) we see that the average over an abelian Chern-Simons field with level 1 repro- 
duces precisely the right sign factors needed to implement the fermionic statistics in 
terms of purely bosonic variables. Inserting this result into the gran canonical parti- 
tion function and using a path integral notation, one can finally prove the bosonization 
identity (2.19). 

When this technique is applied to the t-J model, the electron is described by a 
hard-core boson (/Jq, coupled to an abelian Chern-Simons field. The modulus and phase 
of then interpreted as holon and spinon fields respectively. More precisely, one 

writes 

(^a(z,T) = iy*(z,T)S«(z,T) (2.26) 
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with 

E;(z,t) S«(i,T) = l . (2.27) 
Clearly, there is an arbitrariness in the choice of the phases of H* and Sq,; in fact if 



(2.28) 



the field (pa remains unchanged. However, modulo this U{1) gauge invariance, one 
can still consider H* and Eq, as the generalized "modulus" and "phase" of (fa- Indeed, 
the constraint (2.27) expresses the fact that Sq, is confined on the unit circle like any 
phase, whereas H* must be a hard-core field such that 

H*{i,T) H*{i,T) = 

in order for ip^ to be a hard-core complex bosonic field as required by the bosonization 
procedure. H* and Sq are not yet the most convenient holon and spinon fields, but 
these can be obtained with a further change of variables. In fact, as shown in [18] , there 
are some simplifications if one introduces new fields E{i, r) and Sa{i,T) according to 



-1/2 



E{i,T) 



H*{i,T)^[l-E*{i,T)E{i,T)^ 

S«(z,t)= (l-E*(z,r)E(z,T))"'^' Sa{i,r) 

In terms of these new fields the constraint (2.27) becomes 

^i,T) = E*{i,T)E{i,T) + S*a{i,r)Sa{i,r) -1 = 



(2.29) 



(2.30) 



which is formally similar to (2.11). 

In the abelian bosonization of the t-J model proposed in [17,18], the Chern- 
Simons field A\ is coupled either only to the ho Ions E or only to the spinous 5*0,. 
In the first case the bosonic holons are transmuted into fermions after integrating 
out A\, while the spinous remain bosonic realizing, in this way, the slave fermion 
representation. To see this, let us recall that the bosonized action of the t-J model, 
written in terms E and Sa, is 



S{E,S;A)= rdrl ^ 
•^0 I ren 



E* {i, r) -^Eii, r) + iAo{i, t)[1-E* {i, r) E{i, r) 



d 



+ S*^ (z, r) —Sa (^, r) + iiS*a (z, r) Saih r) 
-t J2 ^*(^,r)e'^<^.>^"^'"'E(j,r)5:(^,r)5«(j,T) + h. 

+ 7ij(^) + 5^iA(i,T)$(i,T)l 



<i,j> 



(2.31) 



10 



where A(i, r) is the Lagrange multipher enforcing the constraint (2.30), and T-Cj{S) 
is given by (2.18b) with ipctii^r) and ^il^l^ii^r) replaced by S^ii^T) and S'*(z,r) re- 
spectively. As mentioned before, in (2.31) the Chcrn-Simons field is coupled (in the 
standard minimal way) only to the holon fields E and E*. Using this action and 
(2.19), the partition function of the t-J model is formally given by 

Q(/5,^) = Q VEVS e-^(^'^;^)^ . (2.32) 

Finally, after computing the expectation value over the Chern-Simons field (or cquiv- 
alently after eliminating A\ through its field equations) one gets the effective action 
of t-J model in the slave fermion representation, i.e. 



E*{i,T)^E{i,T) + S*^{i,T)^So^{i,T) 



Se{i{E,S)= f drl V 

- t J2 [E*{i,T)E{j,T)S*,{i,T)S^{j,T)+h.c.] (2.33) 

<i,j> 

+ HAS) + ^ [/X S*,{i, r) Sail, r) + i A(z, r) $(z, r)] 

Here E and E* denote the hard-core holon fields dressed by the effective Chern-Simons 
contribution. They are defined as follows 

i ^ e{i,j) {E'-{j,T)E{j,T)-i) 
E{i, r) = e E{i, r) (2.34) 

where ©(«, j) is the lattice angle function (see [26,27,28] and Section 4 for a detailed 
discussion of its properties). Using the canonical commutation relations of E and E*, 
it is easy to show that E and E* are fermions satisfying the equal time anticommu- 
tation relations 

[E{i,T),E%j,T)] = S{i,j) . (2.35) 
Furthermore, from (2.34) it is clear that 

E%t,T)E{i,T) = E*{i,T)E{i,T) , 

and hence 

$(i, r) = E*{i, r) E{i, r) + S^h r) S^ii, r) - 1 = $(i, r) . 

Using these results, it becomes evident that the gran canonical partition function 
can be written as a path integral over both fermionic and bosonic fields according to 

Q{P,H)^ J VEVS e-^^"^^^^'^ . (2.36) 
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We observe that this path integral implements correctly the fermionic statistics of the 
original electronic problem. From the effective action (2.33) one can easily derive the 
Hamiltonian of the t-J model in the slave fermion description; when written in the 
operator formalism, such Hamiltonian is exactly the same as the one given in (2.15), 
which we obtained by simply inserting (2.10) into (2.9). 

As remarked in [18], also the slave boson description of the t-J model can be 
realized via abelian bosonization. In fact, if we start from an action similar to (2.31) 
but with the Chern-Simons field A\ coupled only to the spinous Sa, upon elimination 
of Ax, these are transmuted into fermions while the holons remain bosons. Using 
the bosonization formalism, it becomes apparent that the slave fermion and the slave 
boson descriptions are deeply related to each other and can be viewed as two different 
gauge fixings of an abelian gauge invariance exhibited by the bosonized Hamiltonian 
of the t-J model. 

The authors of [18] considered also a non abelian variant of this bosonization 
technique by introducing two Chern-Simons gauge fields: One abelian, and one in the 
fundamental representation of SU{2). In order to represent correctly the fermionic 
statistics of the electronic problem, the coefficients of the two Chern-Simons terms 
must be judiciously chosen (see e.g. Eq. (3.10) of [18]). Furthermore, in order to 
describe only abelian representations of the braid group {i.e. fields whose statistics 
is simply described by a phase and not by a braiding matrix), the level k of the 
SU{2) Chern-Simons term must be related to the spin s of the SU{2) representation 
according to A; = ±2 s. As mentioned, the SU{2) field is taken in the fundamental 
representation with s = 1/2 (the non-relativistic electrons are indeed a doublet of 
SU{2)), and thus k must be either -|-1 or —1. In this non abelian bosonization, the 
abelian field is coupled to the (bosonic) holons while the SU{2) field is coupled to the 
(bosonic) spinous; integrating out these Chern-Simons fields, holons and spinous are 
transmuted into semions. 

Now we are going to generalize (and to simplify) the construction of [18] by 
showing that the semionic statistics of holons and spinous can be realized within the 
context of the abelian bosonization without using non abelian Chern-Simons fields. 
More generally, we will show that holons and spinous in the two-dimensional t-J 
model can be arbitrary anyons subject only to the constraint that the composite 
object made of one holon and one spinon be a fermion. We call this the slave anyon 
representation of the t-J model which, as we will show later, smoothly interpolates 
between the slave fermion and the slave boson descriptions that we mentioned earlier. 
Once again we stress that this "anyonization" is a purely formal operation based on 
mathematical manipulations on the partition function; to see which statistics actually 
occurs, requires a serious investigation of the dynamical aspects of the model. 

Let us now give some details of our construction. We start from the effective 
action (2.33) in the slave fermion picture and, as in [18], we introduce two independent 
Chern-Simons fields, B\ and V\, but, contrary to [18], we take both of them abelian 
^. Then we minimally couple the fermionic holon field E only to Bx and the bosonic 

^ One could equivalently consider the action (2.31) and introduce only a new Chern- 
Simons field V\ coupled to Sa and suitably change the coupling constant in 5'cs(^)- 
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spinon field only to Vx. The resulting action is 



S{E, S; B, V)= jyrl^Y. ^* ^) |:^(^' ^) + ' ^o(^' ^) " ^*(^' ^) ^(^' ^)) 



<»,i> 



,T)e ■'- 



<i,j> 



E{j,T) X 



X 



+ HAS) + S^it, r) S^i, r) + i A(z, r) $(z, r) 



(2.37) 

Notice that the time component Bq is coupled to {1 — E* E) and not simply to the 
density E* E: This is not an artifact but is a direct consequence of expressing the 
original action of the t-J model in terms of E and Sa (cf Eq. (2.31)). The other 
couplings with the gauge fields are no surprise, and one easily sees that 



S{E,S;B,V) 



B=V=0 



= Sefi{E,S) . 



In the spirit of the bosonization procedure, there should be no difference even when 
the two Chern-Simons fields are non vanishing. In particular, the introduction of B\ 
and V\, with their corresponding actions at level ks and ky respectively (see Eq. 
(2.21)), should not spoil the overall fermionic statistics that is already implemented 
in a correct way by the effective action (2.33). Therefore the levels of the two Chern- 
Simons terms must be carefully chosen. Using (2.23), it is not difficult to realize that 
the correct statistics is mantained if 



kB = —kv mod 2 . 



(2.38) 



The minimal choice is 



ks = —kv = — 

V 



(2.39) 



where the parameter v is introduced for later convenience. Contrarily to the non 
abelian bosonization of [18] where the levels of the Chern-Simons terms are both 
fixed, in our approach there is a free parameter which, as we will see in a moment, 
allows to realize anyons of arbitrary statistics and not only semions. 
In this approach the gran canonical partition function is 



Q(/3, /^) = 



(2.40) 
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where the symbol (( denotes the expectation value with respect to the Chern- 

Simons fields B\ and V\ whose levels arc given as in (2.39). Eliminating B\ and V\ 
through their field equations, the resulting action can be entirely expressed in terms 
of new fields, E and Sa and their complex conjugates, defined as follows 

\u {E''{j,T)E{j,T)-l) 

V > ; ^ ' ^ ' 2.41 

-iu ^e(M)s:(j,T)s„(i,T) 

Sa{i,T) = e ^'s" Sa{i,T) . 

These are anyonic fields of statistics v and —v respectively, and, as such, satisfy 
braiding relations among themselves. We postpone the detailed discussion of their 
properties to Section 4; here we simply notice that E is a fermion based anyon whilst 
Sa is boson based. The effective action resulting from (2.40) realizes the slave anyon 
representation of the t-J model and is essentially given by (2.33) with E and Sa 
replaced by E and Sa- However, this representation is purely formal because the 
anyonic fields (2.41) arc non local. In fact, the angle function appearing in the expo- 
nents of (2.41) is multivalued and, to remove all ambiguities, a cut and a fundamental 
domain must be chosen implying that the resulting fields become non-local [29]. Of 
course, if = 1 or = these problems cease to exist since in these circumstances 
the anyons becomes ordinary bosons or fermions. These two cases are clearly in direct 
correspondence with the slave fermion and slave boson representations that we dis- 
cussed earlier, and thus our slave anyon representation smoothly interpolates between 
them. It is worth to point out that the effective action S'eff (£", 3^) is well-defined and 
unambiguous for any value of i/ and does not depend on the detailed definition of 
the angle function; however due to the non locality of E and Sa, it is problematic 
to path-integrate over them. Therefore, in the functional formalism the slave anyon 
representation must be written in terms of bosonic and fermionic fields coupled to 
two abelian Chern-Simons fields as in (2.37), with the gran canonical partition func- 
tion given by (2.40). Things are different in the operator formalism where one can 
define generalized commutation relations for anyonic operators. These are actually 
braiding relations which depend crucially on how the angle function is defined, as we 
will discuss in Section 4. 
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3. The Supersymmetric t-J Model 



When the couphng constants and the chemical potential of the t-J model are 



J = 2t 



II = 7((i) t 



(3.1) 



where 7(d) is the coordination number of the d-dimensional lattice Q, ^ the Hamilto- 
nian (2.1) enjoys very special properties: It can be written as a graded permutation 
operator (up to an irrelevant constant term proportional to the number of lattice 
points) [30] and it commutes with the generators of a superalgebra S'L''(1|2) [6]. In 
d = 1 the t-J model at the supersymmetric point (3.1) can be exactly solved by 
Bethe Ansatz [7,8], and it is tempting to speculate that also in higher dimensions the 
supersymmetry has some non trivial implications. 

Let us first observe that there is a very simple realization of the superalgebra 
SU {1\2) in terms of the projected fermionic operators (2.6) introduced in Section 2. 
In fact, defining 



J+(^)==4(^)c'_(^) , J-(^) = c'l(z) c'+(^) , 
J'{i) = l{n'+{i)-n'_{i)) , 

T(i) = l-i(n'+(i) + n'_(i)) , 
Qa{i) = c'a^) , Ql^) = c'i{i) (a = ±l) , 

and using (2.7), it is immediate to verify that 



(3.2) 



fj+, J-j =2J0 
T,Q 



--Q 
1 



--aJ^ + T 

a 



(3.3) 



Qa 



T,Ql 



1 



where J"*" = X] (''■)■> so on. These are precisely the graded commutation rela- 

tions of the superalgebra S'L''(1|2), which in some circles is also denoted by spl{2,l). 
Notice that Qa and Ql^ for a = ±1 form two doublets of fermionic generators; and 
hence one can say that (3.3) is a N = 2 supersymmetry algebra [31]. 



^ Given any lattice site i E fl, the coordination number is defined as the number 
of points which are nearest neighbors to i. For example, for a square lattice one has 
7(d) = 2d. 
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Choosing the couphng constants of H^j as in (3.1), it is easy to check that 



tJ 







(3.4) 



where is any one of the eight generators of SU {1\2). Thus the t-J model with 
J = 2t and /j, — 'y{d)t is manifestly N = 2 supersymmetric in any dimension d. We 
denote the corresponding Hamiltonian by Hfy^^ . It is interesting to observe that the 
supersymmetry is non-linear ly realized in Hfy^^; in fact the action of the generators 
T"^ on the projected electronic operators yields expressions which, in general, are 
non-linear. For example one finds 



{q+,c'L(z)} = c'L(z)c'+(z) 



(3.5) 



We now show that there is a simpler way to describe this supersymmetry by using 
the slave fermion or the slave boson representations with which it can be realized in 
a linear way. To avoid repetitions, here we limit our discussion to the slave fermion 
case, while in the next section we will discuss in detail the supersymmetry in the slave 
anyon representation. 

Let us then consider the fermionic holon operators e{i) and e^(z), and the bosonic 
spinon operators Sa{i) and sjj,(i) which satisfy the (anti)commutation relations (2.12), 
and let us define 



T = ^(et(i)e(i) + ^4(z).«(z)) , 



Qa = X]e^(i) Sa{i) 



(a = ±1) 



(3.6a) 

(3.6b) 
(3.6c) 



Using (2.12), it is very easy to check that these eight operators close the superalgebra 
SU{1\2) given in (3.3). Furthermore, the action of these generators on holons and 
spinous is always linear. In particular, the SU{2) generators and commute 
with e(z) and e^(z), which are indeed scalars under rotations, while on Sa{i) and sl^{i) 
they act as follows 



f4« 



(3.7a) 



(Of course, in the last two commutators there is no sum over a.) These commutation 
relations simply express the fact that Sq,(z) and sjj,(z) for a = ±1 are two spin 1/2 
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doublets of SU{2). The abelian generator T measures the U{1) charge of holons and 
spinous according to 



T, Sad)] = -^Saii) , [t, 4«] = l^iii) , 

T,e(z)] =-e{i) , [T,et(i)] = e\i) . 



(3.7b) 



Finally, the supersymmetry generators Qa and Ql^ transform the fermionic holons 
into the bosonic spinous and vice-versa. More precisely one finds 



{Qa,eii)} = saii) , {Qi,eHi)} = siii) , 



(3.7c) 



with the other possible (auti)commutators vanishing. 

Before looking at the invariance properties of the Hamiltonian, let us observe 
that the constraint (f){i) in (2.11) is invariant under S'L''(1|2), i.e. 



= . 



(3.8) 



This property is certainly welcome and will have very important implications in the 
following. Let us now turn to the Hamiltonian Htj in the slave fermion representation 
given by (2.15), and chose the coupling constants as in (3.1). A simple calculation 
shows that this Hamiltonian commutes with the generators of SU {1\2) only if the 
Lagrange multiplier A(i) transforms in a non trivial way under the supersymmetries. 
To be precise one finds that 

Qa 1 Htj^ = 



only if 



Qa , A(i)] = -i ^ (a + 7(d)) e\i) Sa{i) 



(3.9) 



where A is the lattice Laplace operator. Even though acceptable on general grounds, 
this transformation property is rather complicated, especially when compared with 
the ones of the other operators entering the Hamiltonian (see (3.7)). Moreover it poses 
serious problems in setting up a mean field theory without breaking the supersymme- 
try. Indeed, the mean field approximation (MFA) instructs us to replace the operator 
A{i) with its mean value X{i) which is a number and, as such, commutes with the su- 
persymmetry generators. Therefore, if A(i) has non trivial transformation properties 
like (3.9), the MFA certainly breaks the supersymmetry. This is very unsatisfactory 
because the mean field approximation is one of the most viable methods to analyze 
the theory in dimensions bigger than one, and if (3.9) were unavoidable, it would 
seem that the role of the supersymmetry of the t-J model can not be investigated in 
this approach. 
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However, this problem can be overcome. Indeed, since the Hilbert space of the 
t-J model is characterized by (f>{i) = 0, we can add to the Hamiltonian (2.15) any 
expression proportional to (f){i) without changing the physical dynamics of the system. 
Thus we can try to modify the Hamiltonian Htj in such a way that the supersymmetry 
be realized without transforming the Lagrange multiplier. This is possible if we 
consider the following operator 



J 



<i,0> 



[i) Sa{i) - e^ii) e(z)) + (^i ^ 



(3.10) 



where Hfj is the Hamiltonian (2.15) with the coupling constants chosen according 
to (3.1). Because of our previous observation, H^j is equivalent to H^j, since their 
difference vanishes on the physical states. With a straightforward calculation, one 
can easily prove that the new Hamiltonian H^j commutes with all generators of 
SU {1\2) without requiring any transformation property of the Lagrange multiplier 
A(i), namely 



and 



A{i) 







(3.11) 



Now the mean field approximation can be used on H^j without breaking the super- 
symmetry. Actually, we can go even further. With simple algebra, the Hamiltonian 
(3.10) can be rewritten as follows ^ 



J 



<i,j> 



+E 



7(rf) (^siii) Saii) + e^ii) e(z)) + i A(z) 4>{i 



(3.12) 



Using (2.11), we then have 
J 



J 



E 4(0^a(i)+e^(0e(j) +$^iA'(^)0(^)--7(ci)l- (3.13) 



<i,3> 



where 



A'(i)=A(i)-i^7(d) , 



and V is the volume of Q, i.e. the number of lattice points. 

This last form of the Hamiltonian is particularly simple and suggestive. In fact, 
using (3.7) it is easy to see that the expression inside the absolute value sign. 



G{hj) = si{i)sc,{j) + e\i)e{j) 



(3.14) 



^ Here and in the following, for any operator O, the symbol \0\ means O. 
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is invariant under SU{1\2). Then, since the constraint (f){i) is also invariant, H^j com- 
mutes with the generators of SU {1\2) only if the Lagrange multiplier A'(z) does not 
transform, as desired. Moreover, when written as in (3.13), the Hamiltonian of the t-J 
model appears as the natural generalization of the Heisenberg XXX model. Indeed, if 
we delete all e and operators from (3.13), we recover precisely the antiferromagnetic 
Heisenberg model which, as well known, is invariant under SU{2). Thus, we could 
say that the supersymmctric t-J model is the natural generalization of the Heisenberg 
model, when the symmetry algebra SU{2) is generalized to SU{1\2). A further gen- 
eralization of the t-J model in this respect, is given by the model recently proposed 
in [31], which is invariant under the superalgebra SU{2\2). In fact, when written in 
the slave fermion representation, the Hamiltonian of [31] is like (3.13) but with two 
sets of fermionic oscillators, corresponding respectively to ho Ions and localons. The 
latter are possible because in this model the double electronic occupancy on the same 
site is allowed and local pairs of electrons (localons) can arise in the spectrum. 

The supersymmetric Hamiltonian (3.13) is particularly suited for a further sim- 
plification: We can perform on it the Hubbard-Stratonovich transformation [32] to 
simplify the J-term, which is quartic in the operators, without breaking the super- 
symmetry. In fact, let us introduce a complex bosonic Hubbard-Stratonovich operator 
Fij defined on the link between the nearest neighbor points i and j, and consider the 
Hamiltonian 



<i,j> 



_ pt p.. 
J ^13 



t J 



+ 



{si{j)s^{i) + e\j)e{i))F,^ 



+ 



J 



^iA'(z)0(z)--7(rf)F 



(3.15) 

Upon elimination of through its "equation of motion" , H'f.j becomes equivalent to 
H'^j and thus it describes the same dynamics. The advantage of writing the Hamil- 
tonian as in (3.15) is that the coupling among the operators is much simpler than in 
(3.13): There are no more quartic interactions. Moreover, since the quantity in (3.14) 
is invariant under SU{1\2) as we have previously remarked, H[j is invariant only if 
the Hubbard-Stratonovich operator does not transform, i.e. 



F- ■ 







(3.16) 



Again this property is welcome, because it allows to make the MFA on H^j without 
breaking supersymmetry. In fact, we can replace the operators Fij and A'(z) with 
their corresponding mean values fij and A'(i), determined in a self-consistent way, 
and get the mean field Hamiltonian 



H 



MF 



= E 



<i,j> 



-jf:jfij+f:j{sl^)saij) + e^{i)eij)) 



+ 



40>a(^) + e'■(J)e(^))/, 



J 



(3.17) 



l{d)V 
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which is invariant under SU {1\2). 

A few remarks arc in order. First of all, once the MFA is made and in particular 
A'(z) is replaced by A'(z), the constraint = is not enforced any more. Then, 
the condition that only the states without double electronic occupancy are permitted 
is released, and in principle new states and excitations, where the number of holons 
and spinous is not restricted as in (2.11), become allowed. Secondly, one can improve 
the MFA without breaking supersymmetry if one considers the phase of the complex 
field fij and the amplitude of X'{i) as fiuctuating. In other words, one writes 

f.. — f . pi Aij 

X'(i) = A(i) + A{{} 

where and X{i) are fixed complex numbers while Aij and A{i) are dynamical 
variables. In particular, if one chooses fij = f and X{i) = A, where / and A are real 
constants, one selects the so-called uniform phase of the slave fermion representation 
[13,20], whose Hamiltonian is ^ 

H = f (4We^^-^a(j) + eUi)e^^-e(j) + h.c. 

+ 5^iA(i) (4(i)s«(i) + et(i)e(i)) 
+ iA J^(4(z)s«(i) + et(i)e(i)) . 

The first two lines can be interpreted as the Hamiltonian of a non relativistic fermion 
(holon) and a pair of non relativistic spin 1/2 bosons (spinous) minimally coupled to 
an Abelian gauge field whose components Ap(i) are 

A^{i)=A{i) , Ai(i) = ^..+i , A^{i) = A,.^^^ (3.20) 

where 1 and 2 denote the unit lattice vectors along the x- and y-axis respectively. This 
gauge field can be identified with the resonating valence bond (RVB) field introduced 
in [2,3,11]. A system similar to that described by (3.19) was also considered and 
analyzed in [33], but there the masses of the fermions and bosons were taken to be 
different, thus breaking explicitly the supersymmetry. 

We conclude this section with a few comments. First of all, the discussion we 
have presented using the operators in the slave fermion representation can be easily 
repeated in the slave boson formalism, where the holons are bosonic and the spinous 
fermionic. In this case most of the formulas, in particular the Hamiltonians (3.13), 
(3.15), (3.17) and (3.19), remain unchanged. The only formal difference is in the 
action of the supersymmetry generators on holons and spinous, since, changing the 

^ We drop from H an irrelevant additive constant. 
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(3.18) 




statistics, the commutators must be replaced by anticommutators and vice-versa. 
Indeed, Eqs. (3.7c) are replaced by 



[Qa, s^ii)^ = e^{i)Sai3 , I^Qi, spii)^ = e{i)Saf3 , 



(3.21) 



On the contrary, (3.7a-b) remain valid also in the slave boson formalism, because J^, 
J° and T are bosonic generators which always close commutators independently of 
the statistics of the operators on which they act. 

Finally, we observe that the analysis of the supersymmetry properties of the t-J 
model as presented in this section, can be equivalently formulated in the functional 
approach, using fields instead of operators. In this case, the action of the generators 
of the superalgebra SU {1\2) on holons and spinous is expressed as transformation 
laws of the corresponding fields, instead of (anti)commutators. For example, using 
the notations of Section 2, the infinitesimal variations of the holon and spinon fields 
turn out to be 



6E{i,T 

SS+{i,T 

SS*_{i,T 



-r]E{i,T) + €+S+{i,T) + €-S-{i,T) , 
+ riE*{i,T) + eXSl{i,T) + elSl{i,T) , 

- ^ 5+ (z, r) - ^+ 5_ (i, r) - I S+ ii,T)-elE{i,T) , 

+ ^ 5_ {i, t) - 9- S+ (i, r) - I 5_ (z, r) - e*_ E{i, r) , 

+ y Slii, r) + e- Sl{i, r) + I 5;(i, r) + e+ E%i, r) , 



y Slit, r) + 9+ Slit, r) + I S^t, r) + e_ i?*(z, r) . 



(3.22) 



Here 9^, 9^ and rj are infinitesimal real parameters associated to the bosonic gen- 
erators J°, J"^ and T respectively, whereas and e* are infinitesimal grassmann 
parameters associated to the fermionic generators Qa and QJ^. It is straightforward 
to verify that the action (2.33) is invariant under the transformations (3.22) if the 
coupling constants are chosen as in (3.1) and the constraint $(z,r) = is enforced. 
Notice that (3.22) are valid both in the slave fermion and in the slave boson rep- 
resentations, the only difference being the statistics of the fields. On the contrary, 
(3.22) are not suited for a discussion of the supersymmetry of the t-J model in the 
slave anyon representation, since in this case holons and spinous are anyonic non-local 
objects for which a functional approach is not well-defined (cf (2.41) and following 
comments). Therefore, the analysis of the invariance of the t-J model under the 
superalgebra 5't/(l|2) in the slave anyon representation must be performed in the 
operator approach, and to this task we turn in the following section. 
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4. The Slave Anyon Representation of the Supersymmetric t-J Model 



In Section 2 we discussed the slave anyon representation of the two-dimensional 
t-J model, and showed that holons and spinons become anyons of opposite statistics if 
one introduces two independent Chern-Simons fields in the bosonization procedure. In 
this section we want to extend the analysis of the supersymmetry properties presented 
in Section 3 to the slave anyon representation of the t-J model. To this aim, we first 
introduce anyonic operators by means of a generalized Jordan- Wigner transformation 
[34,35]. This is given by Eq. (2.41), which, when written in terms of oscillators 
(denoted by lower case letters), becomes 

e(i) = e ^'s" e(i) , 

-i!^ ^ e{i,j) slij) sc(i) 
Sc{i) = e ^es2 _ 

Here e{i) and Sq,(z) are the holon and spinon oscillators in the slave fermion descrip- 
tion, u is the statistical parameter related to the Chern-Simons actions as in (2.39), 
and is the lattice angle function. The exponential factors in (4.1) are some- 

times called disorder operators [35]. One can prove that e{i) and Sa(*) are anyons of 
statistics v and —u respectively, and satisfy braiding relations among themselves. To 
give a rigorous discussion of these braiding relations, a detailed analysis of the lattice 
angle function is necessary. We recall that G(i,j) naively denotes the angle from a 
base point B (eventually taken to infinity) to i as seen by j. However, some care must 
be used to define this angle on a lattice. In fact, to avoid ambiguities, the center of the 
angle 0(i, j) must be a point of the dual lattice near j and a path from B to i must be 
chosen [26,27,28]. Once these two specifications are given, the lattice angle function 
is unambiguously defined. Different choices of the center of or of the path 

along which is measured, lead to different angle functions and hence, through 

(4.1), to different anyonic oscillators characterized by different braiding properties. 
This point was discussed in great detail in [28,36] where the relation between anyons 
and quantum groups was established. Here, instead, we are interested in a different 
perspective, and will choose the angle to define holon and spinon operators in 

such a way that the supersymmetry properties of the t-J model in the slave anyon 
representation become as simple as possible. 

We begin by constructing the Hamiltonian of the t-J model using the anyonic 
oscillators (4.1). We just sketch its derivation without entering the details since only 
elementary manipulations are involved. We start from the action (2.37) and fix the 
coupling constants as in (3.1). Then, we compute the corresponding Hamiltonian, 
rewrite it in the operator formalism and modify it by adding the constraint as we 
did in (3.10). Finally, after removing the Chern-Simons fields, we can introduce the 
anyonic operators (4.1) and obtain the following Hamiltonian 



<i,j> ien 
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where A'(i) is the Lagrange multipher enforcing the constraint </)(z) = to select 
the physical states of the t-J model. This constraint can be written equivalently 
either with the anyonic operators (4.1) or with the oscillators in the slave fermion 
representation; in fact 

0W = 4W Sa{i) + e}{i)e{i) 
= sl(i) Sa(i) + (i) e(i) 

This equality follows from the fact that the disorder operators cancel exactly in the 
combinations sl^{i) Sa{i) and e^{i)e{i), and hence any dependence on the statistical 
parameter i/ and on the angle function drops out from 

For later convenience, and in analogy with (3.14), we introduce the operator 

G{i,j) ^sl(^)sa{j)-e{j)e^(^) for i j , (4.4) 

which appears inside the absolute value of (4.2) and is such that 

\G{iJ)\']=0 . (4.5) 

for all fc, as one can easily check with some straighforward algebra. Notice that 
the ho Ion term of G{i,j) has a relative minus sign with respect to the spinon term, 
and that the creation operator (i) stands on the right of the annihilation operator 
e{j). Obviously the standard normal order can be restored by exchanging e{j) and 
e^(z), but then, due the braiding properties of these oscillators, a phase depending on 
the relative position of i and j is produced (see for example [28]). Of course, when 
u = 0, i.e. when the holons are fermionic and the spinous bosonic, this phase factor 
is simply a minus sign and (4.4) reduces exactly to (3.14). Moreover, in this case the 
Hamiltonian (4.2) becomes equal to that of the slave fermion representation given in 
(3.13), apart from the irrelevant constant proportional to the volume of the lattice 
that we have dropped. 

The anyonic oscillators (4.1) and their braiding properties crucially depend on 
how one chooses the angle function in the disorder operators [28] ; on the contrary the 
Hamiltonian (4.2) is independent of such a choice. To prove this, let us consider the 
operator 

G'{iJ) = S2ii)sM-e'{j)e'\{) (4.6) 

where the primed anyonic oscillators are defined as in (4.1) but with a different angle 
function denoted by Expressing these primed oscillators in terms of the 

unprimed ones, (4.6) becomes 

E [x{i,k)-x{3,k)] 

G'{iJ) = sUi)q'^'' SaU) ^^^^ 



! (4.3) 
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where 
and 



q = e' 



Xihj) 



TT 



(4.8) 



(4.9) 



If we rewrite the exponents in the holon term using (4.3), and then rearrange the 
factors, we get 



G'{i,j) =q 



-X{i,i) at 



4W 



-Yl [x{i,k)-x{j,k)] 4>{k)' 

q-xUd) e{j)e^{i) q 



(4.10) 



From the definition (4.9), it is clear that 

Xih i) = XU, j) = X 



(4.11) 



If we take into account the constraint (f){i) — for all i e fi, the right hand side of 
(4.10) simplifies considerably and one gets 



{ -X+X{j,i)+ [x{i,k)-x{j,k)\ s/3(fc) |- 

G'{i,j) = q^ G{i,j) 



(4.12) 



Therefore, on the physical states of the t-J model, G'{i,j) and G{i,j) simply differ 
by a unitary operator, and thus 



G'{i,3) =G'^{i,j)G'{i,j) = G\i,j)G{i,j)= Gii,j) 



(4.13) 



for all i and j. This equality implies that the Hamiltonian (4.2) does not depend on 
the details of the angle function. This property should not be totally unexpected, 
since different choices of the angle function simply correspond to different gauges for 
the Chern-Simons fields used in the bosonization procedure. 

We now come to the slave anyon representation of the superalgebra SU{1\2). 
As far as the generators J^, J*^ and T are concerned, there is no difference with 
respect to the case discussed in Section 3, and Eqs. (3.6a-b) and (3.7a-b) hold also in 
the anyonic representation, apart from some obvious notational changes. Things are 
different instead for the supersymmetry generators Qa and Ql^. In fact, since these 
are fermionic, they are "sensitive" to the statistics of the operators on which they act. 
For example, in the slave fermion representation Qa and Ql^ close anticommutators 
with e and e ' , but in the slave boson representation they close commutators, as one 
can see from Eqs. (3.7c) and (3.21). Therefore, for intermediate statistics it is natural 
to expect a generalization of these equations. 
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The slave anyon representation of the supercharges can be written in analogy 
with (3.6c); in fact one has 



ien ten 

With simple algebra using (4.1) and (2.12), one can check that Qa and Ql^ in (4.14) 
close the proper anticommutators of SU {1\2), and also that the densities Qa{i) and 
Qlf{i)-i while being fermions with respect to each other, are anyons with respect to 
the constituent holons and spinous. Indeed, Qq(z) and Ql^{i) close anticommutators 
among themselves, but obey braiding relations with e{i) and Sct{i)- Let us now analyze 
if the supercharges (4.14) depend on the choice of the angle function in the anyonic 
oscillators. To this purpose, we consider the density 

Q'aii) = e'Hi) (4.15) 



where the primed oscillators are those introduced earlier. Expressing the latter ones 
in terms of the unprimed oscillators (4.1), we get 

Q'^{i) = e^{i)q s^(i) . (4.16) 

Bringing the Q-factor to the far right and taking into account (4.11), we finally obtain 

Q'aii) - Q"" Qa{i) q . (4.17) 

This equation implies that on the physical states of the t-J model, which satisfy 
= for all i E fl, the two supersymmetry densities Q^(i) and Qa{i) differ only 
by a phase factor independent of i. Thus, on the physical states also the corresponding 
supersymmetry generators simply differ by a constant phase: 

Qa = J2QM = 9^ E^«W =^''Q- , (4.18) 

and hence they close the same superalgebra (3.3). In fact, the (anti)commutation 
relations of SU {1\2) remain unchanged if the supersymmetry charges are rescaled by 
a phase factor. Moreover, it is clear that if Qa commutes with the Hamiltonian, so 
does Q'^. 

The result contained in Eqs. (4.13) and (4.18) is quite important: All possible 
choices of angle function for the anyonic oscillators are equivalent both in the Hamil- 
tonian and in the supersymmetry charges. Therefore, we are free to choose those 
angle functions that are most convenient for our purposes and that simplify most 
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our calculations. Actually, in discussing the supersymmetry properties of the slave 

anyon representation of the t-J model we find very convenient to introduce two angle 
functions that we call horizontal and vertical, and correspondingly define two kinds 
of anyonic oscillators that we also call horizontal (or of type Ti) and vertical (or of 
type V). The reason for this terminology will be evident in a moment. 

For the type 7i anyons, the angle function Qnihj) is defined as follows: First, 
for any point i G O we consider a straight lattice path C^(z) parallel to the x-axis 
from the base point B (taken to infinity in the direction of the positive x-axis) to the 
point i; then we take the dual point 

r = {jx + 1 , + 1) (4.19) 

where and jy are the coordinates of j G O, and £ < 1 is the spacing of a new 
lattice O' in which O is embedded (eventually e can be taken to zero, see [28]). The 
point j* belongs to the lattice dual to O', and can also be viewed as the center of the 
elementary plaquette of Cl' whose lower left corner is j. Finally, we define 0•^(^, j) as 
the angle between B and i measured from j* along the path Cn{i)- This is represented 
in Fig. 1. 

Using this horizontal angle in the Jordan- Wigner transformation (4.1) we get the 
anyons of type 7i, which from now on will be simply denoted by the symbols e{i), 
e'''(z). Sail) and s^{i) with no further specifications. These operators satisfy braiding 
relations among themselves that can be easily obtained using the (anti) commutators 
(2.12) and the properties of the angle function Qnihj)- For example, if i and j are 
two horizontal nearest neighbor points, viz. j — i + 1, one can show that in the limit 



e(z) e(i) + ge(j) e(z) = , (4.20a) 
e{j) e^{i)+qe^{i) e{j) = , (4.20b) 
Sa{i) sp{j) - q-^ sp{j) Sa{i) = , (4.20c) 
^a{j) slii) - q-^ 4(z) Sa{j) = (4.20d) 

where q is given in (4.8). We refer the reader to [28] where these braiding relations 
are discussed and derived in full generality. Here we simply point out that there are 
no phases in the relations between operators in the same point which indeed are 

e{i) e^i) + e\i) e{i) = 1 , (4.21a) 
Sa{i) ~ 4^^) = • (4.21b) 

It is also important to realize that e and Sa commute among themselves; for example 

e{i) sM - Sa{j) e{i) = (4.22) 

for all i and j. From Eqs. (4.20), (4.21) and (4.22), it is straightforward to obtain 
the braiding relations between the supercharge densities Qa{i) and Oa(^) and the 
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anyonic holon and spinon operators. These braiding relations generalize Eqs. (3.7c) 
and (3.21) to the case of arbitrary statistics, and provide an anyonic linear realization 
of the N = 2 supersymmetry. 

Let us now come to the anyons of type V, which are characterized by the vertical 
angle function 6v(i,i). The latter is defined as Qnihj) but with a rotation of 7r/2 
in a counterclockwise way. More precisely, for any point i & Q we consider a straight 
lattice path C\;{i) parallel to the y-axis from the base point B' (taken to infinity in 
the direction of the positive y-axis) to i; then given j e Q, we take the dual point 

*J = {jx - I . + I) (4-23) 

which is the center of the elementary plaquette of Q' whose lower right corner is j, 
and define Q\>{i,j) as the angle between B and i measured from *j along the path 
Cv{i). This is represented in Fig. 2. 

The anyonic oscillators of type V are obtained by using the vertical angle function 
Qvihj) ill (4-1); they satisfy braiding relations among themselves that can be easily 
determined from the oscillator algebra (2.12), and are formally the same as those 
satisfied by the horizontal anyons (indeed, Eqs. (4.20) still hold with j = i + 2). 

These two types of anyonic oscillators are very helpful in simplifying our calcu- 
lations. In fact, if we distinguish between horizontal and vertical lattice bonds and 
correspondingly introduce horizontal and vertical anyons, the Hamiltonian (4.2) can 
be conveniently rewritten as follows 



where 



r I ,2 7, 2 (4.24) 

<i,j>n <i,j>v 



G^{i,j) = sl{i)sM-e{j)e\i) (4.25) 



and G^{i,j) is the same but with vertical anyons. The horizontal Hamiltonian H-j-c 
contains only the interactions between horizontal nearest-neighbor points denoted by 
< i,j >Ui whereas the vertical Hamiltonian ify contains only the contributions from 
vertical nearest-neighbor points denoted by < i,j >v. In some extended sense, H-u 
and Hx> can be considered as one-dimensional Hamiltonians corresponding respec- 
tively to the rows and the columns of the original two-dimensional square lattice. 

Actually, in writing the Hamiltonian (4.24), we could have used only one kind 
of anyonic oscillators, for instance only horizontal ones. But in such a case, the term 
corresponding to vertical nearest neighbor interactions, G'^(z,j), would have been 
more complicated and would have contained extra non-local operators, as shown in 
[27] for a different model. However, as we discussed before, these extra non-local 
operators can be written simply as q raised to suitable combinations of the holon and 
spinon densities, and, after using the constraint = 0, they can be reabsorbed into 
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an overall unitary operator which drops out in taking the absolute value (see (4.12) 
and (4.13)). Finally, we remark that when u = holons and spinous become fermions 
and bosons respectively, and the distinction between horizontal and vertical operators 
ceases to exist. In this case the Hamiltonians H-}^ and i^v in (4.24) can be combined, 
and together they reproduce exactly the Hamiltonian (3.13) of the supersymmetric 
t-J model in the slave fermion representation, up to an irrelevant additive constant 
that we have dropped. 

Let us now discuss the supcrsymnietry properties the t-J model in the slave 
anyon representation. To this purpose, it is convenient to start from the form (4.24) 
of the Hamiltonian and assume that the supersymmetry generators Qa and Ql^ are 
made of anyonic oscillators of type Ti. Since we are free to choose any angle function 
due to the relation (4.18), there is no loss of generality in making this assumption. 
Moreover, to avoid ambiguities we denote the supercharge densities by the symbols 
Q^{i) and Q\^{i)- From the braiding relations (4.20) and the (anti) commutators 
(4.21) and (4.22), it follows that 

Q^{k) , G^{i,j)\ = [Q'^ik) , G^'^ii^j)] =0 ior ij , (4.26) 

and 



+ G'^i^,J) = G-^ihj) (Q^(^) + Q^(j)) 



(4.27) 



Of course the corresponding equations involving the density QU^{i) can be simply 
obtained from (4.26) and (4.27) by hermitean conjugation. These equations imply 
that 

Q«, |G'''(^,J)^J = [Qi, \G^ii,j)\^\=0 , (4.28) 

and hence the horizontal Hamiltonian H-j-i is invariant under supersymmetry. 

To prove that also the vertical Hamiltonian commutes with Qa and Qj^, we first 
observe that Eqs. (4.26) and (4.27) can be translated into the corresponding ones 
for vertical anyons by simply replacing everywhere the superscript Ti with V. Then 
we make use of (4.18) and its hermitean conjugate to rewrite Qa and in terms of 
vertical anyons, and conclude that 



Qa, iG'^iiJ) 



G'^iiJ 







(4.29) 



Finally, it is straightforward to check that the supercharges commute with the con- 
straint operator (/>(i). This completes the proof of the invariance of the whole Hamil- 
tonian (4.24) under the N = 2 supersymmetry. 

In Appendix A we will show that this same result can be obtained by a different 
definition of the angle function in the anyonic oscillators forming the supersymmetry 
generators Qa and Q^^. This alternative method is based on a more intensive use of 
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the constraint (f){i) = 0, but has the advantage of making G{i,j), and not only its 
modulus square, invariant under SU {1\2). 



5. Conclusions 

In this paper we have studied the formal properties of the t-J model in two 
dimensions using the formalism of the slave operators, and showed that, by means of 
a generalized abelian bosonization with two independent Chern-Simons fields, holons 
and spinous may acquire arbitrary complementary statistics in such a way that the 
electron remains a fermion. We would like to stress that the factorization of the 
electron in an antiholon and a spinon is a formal Ansatz, and only a detailed analysis 
at the dynamical level could shed some light on its correctness. 

The t-J model at J = 2t is characterized by the invariance of its Hamiltonian 
under the superalgebra SU{1\2), which is linearly realized on holons and spinous. 
In the case of anyonic statistics, the supercharge densities of SU {1\2) have braid- 
ing properties with the holon and spinon operators, which generalize the standard 
(anti)commutation relations with fermions and bosons. These supersymmetry prop- 
erties of the t-J model at J = 2t could play an important role in the analysis of 
the dynamical aspects of the spin-charge separation, which certainly deserves to be 
investigated further. 
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Appendix A 



In this Appendix we propose two new definitions for the angle function in the 
anyonic oscillators forming the supercharges and that simplify considerably the 
discussion of the supersymmetry properties of the t-J model in the slave anyon repre- 
sentation. In fact with this choice the operator G{i,j), and not only j)| , will be 

supersymmetric. We denote these new angle functions by Qn{hj) ^v{hj) cor- 
responding respectively to the horizontal and vertical anyons. The function Qn{hj) 
is defined as follows: For any point z e f2, we choose the horizontal path Cn{i) from 
the base point B to i, and consider the dual point 

*j = (ix - I ' - I) 

which is the center of the elementary plaquette of O' whose upper right corner is j. 
Then we define Qnihj) as the angle between B and i measured from *j along the 
path Cnii)- This is represented in Fig. 3. 

Using this twiddled horizontal angle in the Jordan- Wigner transformation (4.1) 
we get twiddled anyons of type 7i, which, from now on, we simply denote by the 

symbols e(i), e (i), Sq,(z) and s (i) with no further specifications. These twiddled 
horizontal oscillators are used to define a new anyonic representation of the super- 
symmetry generators. In fact, we can posit 



(A.2) 



As discussed in Section 4, these supercharges are equivalent to those defined in (4.14), 
since they simply differ in the definition of the angle function in the anyonic oscillators 
(cf(4.18)). 

We now study the braiding relations between the supersymmetry generators (A.2) 
and the Hamiltonian (4.24) which is made out of anyons of type 7i and type V. To this 
aim it is therefore necessary to establish the precise connection between Qnihj) ^tnd 
the horizontal and vertical angle functions defined in Section 4. To avoid repetitions, 
we limit our discussion to the horizontal case, but of course similar considerations can 
be made also in the vertical case. 

From the definitions of 0-^(i, j) and Q-j-cih j)i one can derive the following identity 

e„(.i)-e«(i.o={:, £t<i:: 

which holds for any value of e and not only for £ — 0. The derivation of this identity 
is completely similar to that discussed in detail in [27,28], and so we do not reproduce 
it here. Notice that the case i = j is included in this equation since 

Qn{hi) = -^ , G«(i,i) = ^ (A.4) 
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for any i & Q. 

Let us now consider two anyonic holons, one in the horizontal representation and 
one in the twiddled representation. With simple algebra one can prove that 

e(i) ~e{j) = -e-^^-bn{^,J)-enU,^)] g(^) _ (^.Sa) 

Using the identity (A. 3), this equation becomes 

f -Q~^ e{j) e{i) for iy > jy , 
id) ~e{j) = I ^ (A.5b) 
[ -Q e(i) e{i) for iy < jy . 

It is convenient to introduce the function cr, defined on the integers as 

a{m) = 9{m) - 9{-m) + 6^,0 = { \ ^^'^^^ 

where 9{m) is the standard Heaviside step function 

{1 for m > , 

1 for m = , (A.6b) 

for m < . 

Then, we can rewrite (A. 5b) as 

e(i) ~e{j) + q'^^^y-'y) e(i) = , (A.7) 

which is an example of the braiding relations we are looking for. We now consider 
how creation and annihilation operators braid among themselves. For i ^ j we get 



-q e{j) e{i) for iy > jy 
e{i)e'{j)^{ (A.8) 

-g-i e (j) e{i) for iy < jy 



where again we have exploited the identity (A. 3). If z = j instead, using (A. 4) and 
taking into account the inhomogeneous term in the anticommutation relations (2.12a) 
at the same point, we obtain 

e{i) e (i) + q~^ e (i) e{i) = q (A.9) 

where 

x{hj) = -\Qn{hj)-^H{hj)\ ■ (A.IO) 



TT 
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Notice that there is a braiding phase even in the relation at the same point, in contrast 
to (4.21a) which is a standard anticommutator. We can easily combine Eqs. (A. 8) 
and (A. 9) into a single formula as follows 

e{i) e\j) + q-'^^^y-'y) e {j) e{i) = S{i,j) q '^^^ . (A.ll) 

Eqs. (A. 7), (A.ll) and their hermitean conjugates form the complete set of braid- 
ing relations between two anyonic holon operators in the horizontal and twiddled 
representations. 

It is interesting to consider these relations in the limit £ which is part of 
the continuum limit (we recall that e is the spacing of the lattice Q' where Q is 
embedded, and is used to determine the centers from which the angles are measured, 
see for instance (4.2) and (A.l)). The braiding relation (A. 7) remains unchanged 
when £ — 0, whereas (A.ll) simplifies a little, since the function xih k) appearing 
in the right hand side can be computed explicitly. In fact, from the definitions of 
Qn{hj) and Qnihj), it is not difficult to show that, when £ — > 0, 

®nih j) = Onii, j) for iy^jy , 

®nihj) = -T^ , ^nihj)=T^ for iy=jy,ix<jy , (A.12) 

'S>Hihj) = ^niij) = for iy = jy , > . 

Using these values and those in (A. 4), which hold for any e, we get 

x{i,j) = -2Si^jJij,-i,) for £^0 (A.13) 

where 9 is the step function (A. 6b). Then, when £ — > we can rewrite the right hand 
side of (A.12) and get 

e(i) ~e\j) + q-'^^^y-'y) ~e\j) e{i) = 5{i,j) q . (A.14) 

We remark that it is much easier to guess a continuum limit of Eqs. (A. 7) and (A.14) 
than of the analogous braiding relations involving only anyons of one kind [28]. 

What we have done so far with the holons, can be repeated with the spinous. Of 
course, since in the slave anyon representation the spinous are boson based , while the 
holons are fermion based, there are some sign changes; moreover, since the statistics of 
the spinous is opposite to that of the holons, q has to be changed into q~^. Keeping in 
mind these simple modifications, we can immediately translate all braiding relations 
of the holons into those of the spinous. For example, the spinon version of (A. 7) is 

Sail) hU) - q-'^^^y-'y^ ~sp{3) Sail) = . (A.15) 

The other braiding relations of the spinon operators can be obtained in a similar way 
from (A.ll). Moreover it is clear that holons and spinous commute with each other, 
being independent operators (cf (4.22)). 
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After this detailed analisys of the braiding relations, we discuss the supersymme- 
try properties of the Hamiltonian (4.24) using the representation (A. 2) of the super- 
charges. Let us first consider <5a(0 the horizontal nearest neighbor interaction 
G'^ii^j) defined in (4.25) for j = i + 1 (which implies in particular that iy — jy). 
Then, it is easy to prove that 







(A.16) 



for any k ^ The commutation relations of G^{i,j) with the supercharge densities 
in i and j are more subtle. Due to the hard-core condition obeyed by the holon and 
spinon operators in the slave anyon representation, one immediately proves that 

4a(0 , eU) eHi)] = , siii) UJ)] = (A.17) 

On the contrary, after some algebra involving the use of (A. 7) and (A. 11), we get 



, ^ x(i,fc) (et(fc)e(fc)-l) 



^ xU,k) Hk) 



(A.18) 



where in the last line we have transformed the spinon oscillator from the twiddled 
representation to the horizontal one, and have combined all exponents to reproduce 
the constraint operator (f){k). Similarly, one can prove that 



J2 xii,k)Hk) 



(A.19) 



where in the last line the holon oscillator has been transfomed from the twiddled to 
the horizontal representation, and again all exponents have been combined to yield 
the constraint operator. Putting together the last three equations and rearranging 
the right hand sides of (A.18) and (A.19) to bring the factors with the constraint 
operators to the far right, one gets 



ken 

-x{i,j)+ ^ x{i,k) <l>{k) 



(A.20) 



xO',i) + l-xO J)+ ^ xU,k) 4>{k) ' 
q ken 
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From the definition of the function x in (A. 10), it follows that for j = z + 1 



xU, + xU, j) = -xih j) ■ 

Therefore, the numerical factors in the two terms inside the curly brackets of (A. 20) 
are the same and can be factored out, leaving only exponentials of the constraint 
operator. On the physical states of the t-J model where the constraint (fi{i) = is 
identically satisfied, each exponential becomes 1, and thus we can conclude that 



Similarly, one can show that 



Qi,G^{i,j) 







= 



(A.21a) 



(A.21b) 



We now turn to the commutation relations of the supersymmetry generators (A. 2) 
with the vertical nearest neighbor interactions G^{i^j). To simplify the calculations, 
we first notice that all the braiding properties discussed so far, can be translated 
into those appropriate for the vertical anyons by simply using the angles 0v(*,j) 
and Qvihj) in place of Qnihj) and Qnihj) respectively. Notice that is 
defined as Q-uihi) but rotated counterclockwise of 7r/2 (compare Fig.l and Fig. 2). 
Similarly, Ov(^,i) is defined as Qnihj) but rotated of n/2. Then, it is clear that 
all the manipulations we did before naturally apply also to the vertical anyons (both 
with and without the tilde sign). Since on the physical states of the t-J model the 
supercharges constructed by twiddled vertical operators simply differ from those in 
(A. 2) by a global phase factor, we conclude that 







(A.22) 



Finally, since the constraint commutes with the supercharges, from (A. 21) and 
(A.22) we deduce that the whole Hamiltonian of the t-J model in the slave anyon 
representation given in (4.24) is supersymmetric. 

If we compare Eqs. (A. 21) and (A.22) with Eqs. (4.28) and (4.29), we see that us- 
ing the representation (A. 2) for the supercharges, the quadratic interactions G^{i,j) 
and G^{i,j) are themselves supersymmetric, whereas if the supercharges are written 

in the representation (4.14) only the quartic interactions |(5^(z,j)| and |(5^(z,j)| 
are invariant. In both cases, however, it is crucial to use the constraint 4>{i) = which 
selects the physical states of the t-J model. Clearly, in order to check the invariance of 
the Hamiltonian (4.24), the conditions expressed by (A. 21) and (A.22) are sufficient 
but not necessary. However, in view of a possible mean field analysis of the slave 
anyon representation of the t-J model after a Hubbard-Stratonovich transformation, 
it is more convenient to focus on the quadratic interactions G^{i,j) and G^{i,j) 
instead of the quartic Hamiltoinian which derives from them. 
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Figure Captions 



Fig. 1 Representation of the lattice angle O-j-cihj) between the base point B and i 
measured along the horizontal curve Cuii) from the point j* . 

Fig. 2 Representation of the lattice angle between the base point B' and i 

measured along the vertical curve C^{i) from the point *j. 

Fig. 3 Representation of the lattice angle Q-niiij) between the base point B and i 
measured along the curve C-uii) from the point *j. 
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Abstract 



We discuss the properties of the supersymmetric t-,J model in the formahsm of the 
slave operators. In particular we introduce a generalized abelian bosonization for the 
model in two dimensions, and show that holons and spinous can be anyons of arbitrary 
complementary statistics (slave anyon representation). The braiding properties of 
these anyonic operators are thoroughly analyzed, and are used to provide an explicit 
linear realization of the superalgebra SU{1\2). Finally, we prove that the Hamiltonian 
of the t-J model in the slave anyon representation is invariant under SU {1\2) for 
J = 2t. 
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